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Abstract This paper deals with the validation and theoretical analyses of a simple method for 
calculating the intrinsic rate of increase, rm, of aphids and mites. This method does not require a 
detailed fecundity table data. The value of rm can be estimated by the simple equation: rm=0.74 ln 
(Md)/d, where d is the time from birth to first reproductions Md is the number of female offspring 
produced per original female in 2d. This method was developed by Wyatt and White. We reached the 
following conclusions from our study: When the parameter, Md is less than l, this equation is not 
appropriate for estimating rm of populations, When the parameter, Md is larger than 1 and about 
70% or more of the reproductive contribution to the rm is achieved in 2d, the equation is appropriate 
for calculatng the rm of aphids and mites.As Md decreases, the required reproductive contribution 1o 
the rm in 2d corresponding to the constant 0.74 will increase. However, whether 70% or more of 
reproductive contribution to rm is achieved in 2d still is a problem for whole taxa of aphids and 
mites. Therefore, further study on reproductive distribution will be required for practical application 
of this method, 
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introduction 


The intrinsic rate of natural increase, rm, is one of the most important biological properties of asimal 
populations. Birch (1948) developed an accurate method of calculating rm on the basis of Lotka's integra] 
equation (Lotka, 1925). Thea, approximate estimate methods were developed with a view to simplifying 
calculation of rm values (Leslie ef al., 1940; Birch, 1948; Howe, 1853; Laughlin, 1965. However, 
these methods depend on detailed fecundity table data. In many cases, a fecundity table is very difficult 
to compile. Fortunately, much attention had been paid to the relationships between the rate and more 
readily obtained parameters (Cole, 1954) Morisita, 1965, Fenchel, 1974) Wyatt ef al., 1977) Parise, 
1978; Gaston, 1988). These studies have resulted in some simple methods io estimate the intrinsic rate of 
increase. However, the wide application of these methods depends on how they work in practice. 
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Wyatt & White (1977) developed a simple method to determine the rm values of aphids and mites 
based on their reproductive distribution. This method simplifies the calculation procedure, ard on the 
other hand, does not require the detailed fecundity table. In their esse, this equation worked very well 
for nine species of aphids and four species of tetranycbid mites. Gerson ef al (1883) also obtained close 
agreement between the rm calculated by Wyatt & White method and the accurate rm calculated by the 
Birch method for the acrid mite Rhizoghyphus robini. In this paper, we cea] witb the fur her validation 
of this metbod, and furthermore, discuss the application conditions under whieh i. cen provide a god 


estimate of the rm value. 


Materials and Methods 


Fecundity tables of the winter populciion srd summer populelien of woolly apple apbid, Eriosoma 
lanigerum (Hausm} were constructed fram leboaralory observations cerried out in Kunming, Yunnan 
province of China. Apple twigs were collected from Nanlongtan Apple Orchard in the Kunming area in 
January and July of 1987, respectively. Aphids were resred on these twigs under controlled conditions- 
Each laboratory population consisted of 600—890 apbids distributed over 40 twigs which were placed in 
8 beakers on the average. Life history observations and records were made daily from the birth of apbids- 
The other fecundity tables used in tbis paper were extracted from published literature. 

The intrinsic rates of natural increase were calculated by the Birch method and the Wyatt & White 
method. The following iwo equations were developed by Birch (1948) to calculate rer on ihe basis of 


detailed fecundity table data. 


> exp (-1rmx)? leamz= 1 cena CJ) 
Zod 
ta=in(Ry)/T mhmmspmenens (2) 


Where x is the ages la and mæ are age-specific survival rate and age-specific fecundity, respectively 
Ro and T are the net reproductive rate and meen generation time, respectively, Ro and T are given by the 


following ¢epuations: 
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Equation | is the accurate method for calculating the rale rug Equation 2 is an approximate method fur 

estimating the rate rm. The equation developed by Wyatt & While (1977) takes the following form: 
tm=c In(Md) /d ama CGG 

where d is the prerepreductive time from the birth ‘o the.first reproductions Md is the number of female 

offspring produced by per original female in 2d, measured from birth. When the first day of birth is 

designated as zero, Md could be expressed as the following: . 


aS Zd-—t 
Md= 2) lomz= >) ix/No forle=Na/No and mz=fa/Ne eneee (6) 


where NO is the number of original Iemales; Ne is (he number of females alive at the age of xa Ia is 
the number of female offspring at the age of x, In equation 5, the constant c is a correction factor- 
Wyatt snd White (1977) found out ihat ihe value of c was around 0.74 on the basis of regression of 
accurate om on the factor In¢€Md} /d. 

tor= 0.74 In€Md)/d ee D 


Rearranging equation | will produce the cquation &. 
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$l exp (-rmx) lame+ $] exp C- rmx) lemz= I T sementa ABI 
m@=d az=23d 


According lo the method developed by Birch (1984), the value of first item on the left of ooa on & can 
be dcfined as the reproductive contribution to the rm during 2d. i O 

In addition, the model without a constant item, shaped like y =ax, is used in the following regression 
analyses of the accurate rm on the approximate estimates of rm. The AST 386/16 compuler and software 
SYSTAT 4.4 were used in this study. 


Results and Dicussions 


1. Comparison of Birch and Wyatt & White methods 

Since both the accurate and approximate Birch methods are widely used in calculating rm of popula- 
tiors, it is necessary to compare a new method with the Birch methods. Substituting the fecundity table 
data into equation 1 and 2, we could calculate the accurate rm values and the approximate rm values 
(Table 1). The regression analysis suggested that the epproximate Birch method could provide a good 
estimate of the rm value with a standard error of estimate of 0.021 (Fig. 1). Correlation coefficient of 
the approximate rm values with relation to the accurate rm values was 0.996 (P<0.001). 

The parameter, Md, was calculated by equation 6 (Table 2). The regression of the accurate Birch 
rm values on the factor In (Md)/d produced the equation 9. That is, 

tm= 0.691 In (Md) /d rerecrearsee {BY 

The correlation coefficient between them was 0.938 (P<0.001). This implies that the linear correlation 
is acceptable. However, estimate of rm by equation 9 produced greater value of the standard estimate error 
compared “with the approximate Birch method estimate (Fig. 2). The fomer was 0.075, while the latter 
was 0.021, Similarly, we directly used the equation, rm=0.74 ln (Md)}/d, (equation 7), proposed by 
Wyatt and White, to estimate the rm. The resulting estimated values of rm still had a standard error of 
estimate of 0.079 from the accurate rm (Fig. 3). 


Tab. 1 intrinsic rates of increase Crm) computed by Birch methods 


TS 





Species Experimental conditions on cm Source 
Eriosoma lanigerum (Srring) 20% 7o%RH® 0.1644 0.1390 Kuang at al., 1989 
E. lonigerum (Summer) 20% 70%RH 0.1675 0.1642 Original 
E. lonigerum (Autumn) 26% 70%RH 0.2952 0.2848 Shan et al., 1988 
E. lanigerum (Winter) 200 70%RH —0.0565 -0.0591 Original 
Eotetranychus hicorice 21.10 24:0 LDG 0.0771 0.0750 Micinski ef af., 1981 
E. hicoriae 21.10 24:0 LD -0.0001 -0.0001 Same 
E. hicoriae 26.00 24:0 LD 0.0974 0.0928 Same 
E. hiicorae 32.2% 24:0 LD 0.1051 0.1027 Same 
Myzus persicae — 9.4209 0.3227 Harrisor, 1969 
Rhopalosiphum maidis —— 0.4041 0.3081 Wu ef al., 1986 
Zetzellia mail 90 56%RH Q. 1084 0.1072 White & Laing, 1977 


Notes: @ The rm values estimated by the accurate Birch methods © The rm values estimated by the 
approximate Birch methods @RH is relative humidity, © The ratio of light to darkness is 142 (0 for imm- 
aturre stages. 


ni ma mhana ee te re. ot eo 


OOO Obe SASS SE A dE I A a a mA me PTY * 


ar nears aoa 


ào awe se RS ies 





Tab. 2 Th values of parameters associated with the simple method 


Species Experimental § Prereproductive Md Reproductive contribution 





conditions time, d (days) to the rm in 2d, Cd 

Eriosoma lanigerum (Spring) 200 70%RH® g 6, 2465 0.800: 
E. lonigerum (Summer) 200 70%RH 5 0.8352 0, 2028 
E. fanigerum (Autumn) 260 70%RH 7 19.4752 0.8250 
E. lonigerum (Winter) 200 7o RH 6 0.1772 0.2918 
Eotetranychus hicoriae 21.10 24:0 LDG 25 13.3629 0.9995 
E. hicoriae 21.10 24:0 LD 15 0.8208 0,823! 
E. hicoriae 26,00 24:0 LD iL 4.1910 0.8919 
E. hicoriae 32.20 24:0 LD ll 4.4420 0,9961 
M yzus persicae — 5 13.8740 0.9277 
Rhapalosiphum maidis — 4 9, 2820 0.6824 
Zetzellia mali 190 56CRH 16 10,0407 1 





Notes: @ RH represents relative humidity. @ The ratio of light to darkness was 14:0 LD for immature 
slages. 


In order to explain these results, we calculated the reproductive contribution, Cd, to the rm in a period 
equal to 2d from the birth by substituting the accurate Birch rm and fecundity table data into the first 
item of equation 8. 

Equation 8 clearly suggests that reproduction after 2d can be ignored only when the first item in this 
equation accounts for sufficiently high proportion. In our instance, the reproductive contribution to the tm 
in 2¢ were so low for summer population and winter population of E. fonigerum that the reproduction 
after 2d could not be ignored (Table 2). Therefore, application of this method is not appropriate in this 
situation. 

2. Theoretical analyses of Myatt & White method 

The firs} item on the left side of equation 8 represents the reproductive contribution to the rm in 2d. 
When its value equals ], that is, reproduction after 2d does not occur, equation 8 is simplified into 


equation l0, 
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Tn by the approximate Birch method 


rig. 1 The regression of the rm values computed by the accurate Birch method 


on the rim values computed by the approximate Birch method 
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Fig. 2 The regression of the rm values computed by the accurate Birch met hod 
on the rm values estimated by the equation; rm=9.601 In (Md)/d 
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Ta by the equation: ra = 0.74 log(M,y/d 
Fig. 3 The regression of the rm values computed by the accurate Birch method 
on the rm values estimated by the equation; rm=0.74 in (Mdy/d 
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In this case, the accurate rm can be calculated according to this equation, while the approximate rm can 
be estimated according to the approximate equation developed by Birch. That is, 
Resexp CrmT) or exp C- rm Rosi tts (41) 


Where T is generation time. in equation II, the net reproduction rate, Ro, can be calculated by the 
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following equation: 


S ~~. = aterm pmnan o a as es rea e A a a a 
Ro = lems 
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Apparently, Ro in this case is equal to the number of female offspring produced by per original female in 

2d, Mg, according to equation 6. That is, 
Ry = Md 

Substituting Md for Ro in equation H. produces the following cquation; 
Md = exp CrmT ? or exp C- re Th = l a sernenaessa {F2 

Comparing equation 30 with equation 12, the Tollowing equalion will saturally resull, S h 
cxp (-rxT) Mg we St exp — mx) lume. <-1 ages pine: ues Samah tases C14) 
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When reproductipn after 2d occurs, we should add the cond: item of oust 8 io the two sides of 
equation 13. Tat is, ` 


wW 
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Where Ca is the reproductive contribution to the rm in 2d (Table 2}. As a result, tm is expressed as 


follows: 7 
rm=in (Ma/Cay/T ~~. 
Le! T =d/Co, where Cy is the m thuss : 
rm=Co in (Ma/Ca)/d ag s aie * (15) 


Comparing equation 15 with the Wyatt & White “equation equation 5), it is discovered that when the 
reproductive contribution to the rm in 2d, Cg, is equal to Í or near } (e. gk. 9.95 or 95%), equation 15 
becomes equation- 5. Apparently, equation 5 is a special. case of this.equation.. 

We use equation 15 to analyze data in the case of this study. The- regression of i accurate rm by 
Birch on the factor [n-<Mz/Cg)/d-gave regression- soefficient-of- 0.676.. abst aes 

te= 0.678 In (Ma/Ca)/d - DU wanasens ons (16) 

The correlation coefficient was 0.995 (P>0.001), “The estimate of rm by equation 16 had a standard 
estimate error. of 0.023., Compared with. SRS: 9s k estimated., iiai af. rm wera more accurate 
(Fig. 4). pos 45 

Unfortunately, we cannot obtain the siimaa alee of Cai in the Daed detailed fecundity table 
data or the accurate rm. Therefore, equation !5 or 16 cannot replace equation 5 developed by Wyatt. & White. 
However, Equation 15 can Provide some important information on the simple method developed: ‘by Wyatt 
and White and further development of this method. Rearranging equation [6 will — the following 


cquation: 


rm= {Co In (Md/Ca}/In May} log Madd Ss Shanen (1P) 
Comparing equation {7 with equation 5, we deduce ihe following equation: . C 
C=C, In (MaC ina) l terrearen iaa C18) 
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"Since Cp = 0.876 in equation l6, 
Thus, C=0.676 In (Ma/Cad/in (Md tt ttre C199 
Equation 18 and 19 are two important equations. Here, they can be used to analyze how the constant c 


in equation 5 changes. with. wearying- parameters Ca and Mags- and to. seek forthe conditions under which 
the constant c will have” a valne around 0.74. 





According to cquation 19, results of plotting the- constant e Po the reproductive contribution to 
the rm in 2d, (Ca), suggest that the constant c is around 0.74 for Ma>1 and Ca>0.7 Cor 70%) 
(Fig. 5). In the case of our study, 8 sets of values of parameter, Ma, were larger than l and ranging 
from 4 to 20, and the reproductive contributions to the rm in 2d, Cd ranged from 0.6824 to ] (68.24% 
to 100%). Therefore, we theoretically predict from Fig. 5a that the c in equation 5 should be around 
0.74 for these $ sets of data. The regression of the accurate rm by Birch on thp factor In “Md)/d frem 
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only these eight sets of data, suggested that the constant e was 0,747, That is, 
tm= 0.747 In (Md) /d teeeratesere C20) 
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Fig. 4 The regression of the rm vaites “computed by the accurate Birch method 
on therm valties estimated by the equation; .rm= 0.876 In (Md/Cdy/d 

The correlation Sativa: -was -as high as 0.999 (P<0.00)). The standard error of estimate was only 
0.014. This is in complete agreement with our prediction, ard shows that estimates “of the tm calculated 
by <quation 20 are more accurate than apprceximate Birch method estimates (Fig. & 63. This suggests 
that, 9.74 is an appropriate value of constant c in equation 6 for Cd>0.7 Cor 70%) and Md>1. In addition, 
the value of Cd in 2d corresponding to 0.74 will increase, as the parameter Md, decreases (Fig. 6a). 

When the parameter, Md, is less than 1, the constant 9,74 is almost not ‘Appropriate unless the repro- 
ductive contribution to the rm in ad is s virtually 1 ce. 8- 98%, n até } vies 5b). When eee 
of populations. According.to, the equation 19, values of the constant c were 0.1948 and 0.0098 for the 
winter population of E. lenigerum and. ore. population of E. hieorige (21.10, 24:0LD) » respectively. 
In addition, if the reproductive contribution to the tm in 2d is too low, the Md can be less than } for 
positive” growth of populations. In this case, 9.74 is not appropriate, either. For example, the constant 
g for summer population E. lanigerun, haying an Md of 0.8152, was ~ 4,6024. 
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Fig. 5a The effect of the parameter Md and reproductive contribution to the rm in 
2d on the correction factor c Le =0.676 In (Md/Cd)/lIn(Md)) 
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Fig. 5b The effect of the parameter Md and reproductive contribution to the tm in 
24 on the correction factor ¢ C¢ =0,676 In (Md/Cd}/In (Md) J 
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Pa by the equation: ra = 6.747 log(M,)t 


Fig. 6 The regression of the rm values computed by the accurate Birch method on 
the rm values estimated by the equation; rm=0.%47 In (Md)/d 


Conclusions 


When the parameter, Md, is less than l, this equation, rm=0.74 In Md/d, is not appropriate for 
estimating rm of populations. The Md of less than | may result from too low reproductive contribution 
to the rm in 2d or negative growth of populations. [n these cases. the cons tant c does not closely correspond 
to 0.74. Of course, when the reproduction contribution to the rm in 2d can be considered to be ] œ. g. 
98%, 99%), this equation still can provide an approximate value of rm. 

When the parameter, Md ¿the number of female offspring produced by per original female in 2d from 
birth, where d is the time from birth to first reproduction) is larger than 1 and 70% or more of contri- 
bution to rm is achieved in 2d, the equation developed by Wyatt and White (1977), is appropriate for 
calculating rm of aphidsTand mites. As Md decreases, the required reproductive contribution toZthe tm in 
2d correspording to 0.74 will increase. 

However, it must be stressed that we could oot determine the reproductive contribution to rm in 2d 

without detailed fecunditty table data. Delosch (1974) found out that about 0.86 Cor 96%) of reproductive 
contribution of aphids was achieved in 2d. But whether this is true for the whole taxa of aphids and the 
other insects remains to be further studied. This is quite important to practical application of Wyatt-Whiie 
method. 
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BARRE, HA rm WARM: rm=0.74 In Mdd. OH, d SRE, WANEN- KEAS J 
ANRE, Mok wR Ads ST EER. KOR Watt A Whitedt 09, RER BT 
SRA: SMdht 1, RR WRRARA, HMAF 1 Medd ee REM AS IO“, TURR 
AHEHE rms 2 Rak, MA RAA HREM rm RM 0% RTO LLB 
PAA. TESLA AH SEA Fa A Ha A 


